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Abstract 

We introduce a new method for determining the critical indices of 
the deconfinement transition in gauge theories. The method is based 
on the finite size scaling behavior of the expectation value of simple 
^ ■ lattice operators, such as the plaquette. We test the method for the 

case of SU(3) pure gauge theory in (2 + 1) dimensions and obtain a 
precise determination of the critical index v, in agreement with the 
prediction of the Svetitsky-Yaffe conjecture. 
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Gauge theories at finite temperature undergo a phase transition from a low 
temperature phase, in which color sources are confined by a linearly rising po- 
tential, to a high temperature phase in which the attractive force is screened 
at long distances and isolated quarks can exist. This same behavior can be 
observed even in pure gauge theories, without matter fields, by studying the 
theory in presence of static sources. 

The order parameter of the deconfinement transition in a pure gauge 
theory is the Polyakov loop [|], |2|, representing the world-line of a static 
source of color charge (that can be thought of as an infinitely massive quark) . 
The Polyakov loop takes values in the center of the gauge group: the (d+ 1)- 
dimensional gauge theory generates an effective <i-dimensional theory for the 
Polyakov loop, having the center of the gauge group as global symmetry 
group. The deconfinement phase transition coincides with the spontaneous 
breaking of this global symmetry. 

In this work we introduce a new method for the computation of the 
critical indices of a pure gauge theory with second order deconfinement tran- 
sition. According to the Svetitsky-Yaffe conjecture |J, such critical indices 
for the (d+ 1) -dimensional gauge theory coincide with those of the effective 
(i-dimensional model, if the latter has also a second order phase transition. 
We will consider in particular the critical index v of the correlation length. 
For a general d-dimensional statistical model a possible way of extracting 
the value of v from lattice Monte Carlo simulations is to study the finite size 
scaling (FSS) behavior of the energy operator: FSS theory predicts that, if 
L is the lattice size, the expectation value of the energy operator behaves for 
large L as 

(E)l ~ (E)oo + kL»~ d , (1) 

where (E)^ is the expectation value of the energy operator in the thermo- 
dynamic limit and A; is a non- universal constant. This method was applied 
e.g. in Ref. to evaluate the critical index v of the three-dimensional Ising 
model. 

Compared to other methods based on the FSS of fluctuation operators 
such as susceptibilities or Binder cumulants, the advantage lies in the fact 
that (E)l can be computed to high accuracy. The main drawback is that 
the term containing v in Eq.([I]) is subdominant for L — > oo with respect to 
the bulk contribution (E)^. The numerical results we will present strongly 
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suggest that, in the case of gauge theories, the advantages outweigh the 
drawbacks and the method can give very accurate results. 

Consider a gauge invariant operator O that is invariant also under the 
global symmetry given by the center of the gauge group. Examples of such 
operators include Wilson loops and those of the form 

P{x)P\y) , (2) 

where P(x) and P{y) are Polyakov loops at two different sites x and y of the 
(i-dimensional space. It is natural to expect the operator product expansion 
(OPE) of any such operator to have the same form as the one of E: 

O = Cl I + c e e + - ■ ■ . (3) 

Here / and e are, respectively, the identity and the scaling energy operator[| 
in the statistical model, and the dots represent contributions of irrelevant 
operators. The ansatz of Eq. (|3|) was introduced and tested in Ref. ||, and 
used in Ref. || 0], to obtain some exact results on correlation functions of 
(2 + l)-dimensional gauge theories at the deconfmement transition. 

In particular, Eq. (^) implies that the FSS behavior of the expectation 
value (O) will have the form of Eq. (|]): 

{6) L ^ {0) oo + cL 1 -- d . (4) 

The contributions of the irrelevant operators will be subleading for L —>■ oo. 

We conclude that the FSS behavior of any such operator can be used to 
determine the value of v through Eq. (f4|). The obvious advantage is that 
one can use operators, such as the plaquette, whose expectation value can 
be computed to high accuracy with relatively modest computational effort. 

To test the method in practice, we used it to determine the critical index 
v for SU(3) pure gauge theory in (2 + l)-dimensions. This was evaluated 
from Monte Carlo simulations in Ref. |§ with the result 

v_uc = 0.90(20) . (5) 

-'^Not to be confused with the lattice energy operator E, which is just an example of an 
operator with OPE given by Eq. (fj). 
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However, if the Svetitsky-Yaffe conjecture || holds then the 2-dimensional 
effective theory of the Polyakov loop is in the universality class of the 3-state 
Potts model, and the value of v is known exactly: 

u=- = 0.833- •• . (6) 
6 

The observables we considered are the expectation values of the plaquette 
and of the operator defined in Eq. @, where x and y are taken to be nearest 
neighbors sites in the 2-dimensional (spatial) lattice. We computed these 
expectation values on lattices with temporal extension N t = 2 and spatial 
sizes L = N x = N y ranging from 7 to 30. The simulations were performed at 
the critical coupling 

(3 c (N t = 2) = 8.155 (7) 

obtained in Ref. H. The simulation algorithm we adopted was a mixture of 
one sweep of a 10-hit Metropolis and four sweeps of over-relaxation consisting 
in two updates of (random) SU(2) subgroups. For each simulation we col- 
lected 400K equilibrium configurations, separated each other by 10 updating 
steps. The error analysis was performed by the jackknife method applied to 
data bins at different levels of blocking. 

We report in Table 1 the expectation values we obtained. Because of the 
asymmetry of the lattice, space-like and time-like plaquettes have obviously 
different expectation values and must be considered as two different oper- 
ators. To evaluate u, we performed a single multibranched fit of the three 
data sets. To avoid cross-correlations we included in the fit only space-like 
("magnetic") plaquettes from lattices with odd L and time- like ("electric") 
plaquettes from lattices with even L. Polyakov loop correlations were mea- 
sured in separate runs and therefore are not correlated with the plaquette 
measurements. 

The result of the fit is 

v = 0.827(22) , = 0-84 , (8) 

in excellent agreement with the value of Eq. (||), coming from the Svetitsky- 
Yaffe conjecture, and with remarkably improved accuracy in comparison to 
the existing Monte Carlo evaluation. 

In conclusion, the method we have proposed appears to give very precise 
results for critical indices, thanks to the fact that only "simple" expectation 
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values have to be evaluated, such as the plaquette expectation value. The 
same method can be used to evaluate the critical coupling (that, in our case, 
was taken from the literature): since Eq. (^) is valid only at the critical point, 
one should perform the same fit at several couplings and determine the value 
of the critical one, by comparing the x 2 's of the different fits. 
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Table 1: Expectation values of the electric and magnetic plaquette and 
Polyakov loop correlator. 



L 


Plaquette (el.) 


Plaquette ( 


mag.) 


PF correl. 


7 


0.628855(23) 


0.625721 


'22) 




8 


0.628562(30) 


0.625494 


'27) 

\ } 


1.2707(13) 


9 


0.628318(26) 


0.625358 


'24) 


1.2530(15) 


10 


0.628107(27) 


0.625180 


'25) 


1.2358(15) 


11 


0.627929(26) 


0.625059 


'21) 




12 


627825(25) 


624983 

V 7 • V7 -L '.7 v 7 '.J 


'21) 


1.2122(14) 


13 


0.627729(26) 


0.624893 


'22) 


1.1997(15) 


14 


0.627593(23) 


0.624821 


;2o) 


1.1944(16) 


15 


0.627493(24) 


0.624741 


19) 




16 


0.627416(22) 


0.624692 


19) 


1.1829(15) 


17 


0.627343(23) 


0.624643 


19) 




18 


0.627270(22) 


0.624598 


!l7) 


1.1669(16) 


19 


0.627212(23) 


0.624545 


18) 


1.1625(16) 


20 


0.627182(22) 


0.624523 


:i7) 


1.1573(16) 


21 


0.627121(23) 


0.624482 


19) 




22 


0.627067(20) 


0.624445 


16) 


1.1511(15) 


23 


0.627007(23) 


0.624396 


:i7) 




24 


0.627043(23) 


0.624426 


;i7) 


1.1429(16) 


25 


0.626977(23) 


0.624385 


;i7) 


1.1408(17) 


26 


0.626922(22) 


0.624338 


16) 


1.1390(16) 


28 


0.626855(24) 


0.624297 


;i7) 


1.1322(16) 


29 


0.626869(22) 


0.624311 


;i7) 


1.1311(19) 


30 


0.626814(26) 


0.624265 


18) 


1.1289(17) 
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